Formation of Vector Solitary waves with Mixed Dispersion in Bose-Einstein 

Condensates. 
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We demonstrate the existence of a new class of two-component vector solitary waves in which 
dispersion coefficients have of opposite signs. Stability is achieved by inclusion of an additional 
linear coupling between the vector components that counterbalances the instability produced by 
the mixed dispersion and the non-linearity. In addition, we demonstrate that these solutions are 
experimentally observable as gap vector solitons in Bose-Einstein condensates located in oscillating 
optical lattices. The proposed experiment falls well within the standard experimental procedures 
for generating gap solitons in today's laboratories. 
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I. INTRODUCTION 

Dispersion is a fundamental process in the dynamics 
of wavepackets. Its effect is revealed dramatically as the 
expansion of the wavepacket's size in time. For practical 
applications this extension becomes a serious limitation: 
degrades communication bit rates, and prevents efficient 
methods for storage and transport of classical or quantum 
information. Fortunately, this situation changes radically 
in non- linear systems, where the non- linear potential can 
counteract dispersion and may result in a preservation 
of the wavepacket shape in time. The characterization 
of situations where these soliton solutions emerge consti- 
tutes a fundamental challenge for research. Besides its 
mathematical interest, especially relevant is the descrip- 
tion of situations in which solitons or soliton-like struc- 
tures appear in Nature, opening the door for technologi- 
cal applications (as happened with electromagnetic waves 
in Kerr media 

The non-linear Schrodinger equation (NLSE) consti- 
tutes one paradigmatic case for the study of solitons 
0, 0. It describes different physical systems, among 
them the mean field dynamics of Bose-Einstein conden- 
sates (BEC). Bright solitons, or more generally solitary 
waves, in BEC have been demonstrated experimentally 
for the case of attractive interactions 0, || and dark soli- 
tons in the case of repulsive interactions 0, • In ad- 
dition, bright solitons can be also obtained for repulsive 
interactions in condensates near the edge of the Brillouin 
zone of an optical lattice (gap solitons) . In this case, 
a condensate in the lowest band will experience a nega- 
tive effective mass, thus inverting the sign of dispersion, 
which now can be counterbalanced by the repulsive non- 
linear potential. The first experimental observation of 
this type of solitary waves has been reported recently @ . 

Vector solitons may appear in situations where two or 
more different wave modes can coexist. In this case, the 
soliton is composed by a set of wavepackets, one for each 
mode that propagate jointly. Every mode, therefore, con- 
tributes to the overall non-linear potential that allows for 
self-trapping. Although the integrability of the resulting 



system of equations is not ensured for the general two- 
component case, special cases as biased photorefractive 
crystals can be demonstrated to enclose Manakov soli- 
tons as solutions Other examples of two-component 
vector solitary waves can be found in waveguide arrays, 
where the two radiation modes correspond to different en- 
ergy bands flol ITT| . In the general case, each equation of 
the set describing the vector soliton includes the counter- 
balance between dispersion and non-linearity. Therefore 
each vector component, even evolving uncoupled, has 
the dynamical conditions to reach a soliton wave, which 
grants in great manner the stability of the coupled sys- 
tem. However, this condition is quite restrictive, since 
one may think in situations in which not all the vector 
components are described by equations in which the non- 
linearity and the dispersion balances. A particular case 
is the formation of vector gap solitons in which the geom- 
etry of the bands at the edge of the Brillouin zone only 
leads to a negative sign in the dispersion of the lowest 
band. In the uncoupled case, therefore, only the vector 
component corresponding to the lowest band will evolve 
to a soliton, while the other components (belonging to 
the upper bands) will rapidly disperse. In this case, the 
coupling is a fundamental aspect for the formation of a 
localized vector wavepacket. The aim of this letter is, 
first, to propose for the first time a set of two coupled 
NSE with opposite signs in the dispersion term which 
sustains localized self-trapped vector solutions, and, sec- 
ond, to demonstrate that such solitary solutions can be 
observed in physical systems as BEC in oscillating optical 
lattices. 

To begin, let us consider the following system of cou- 
pled NLSE equations: 
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Where a, 71 and 72 are positive, a accounts for the 
mass difference between the two vector components, Aft 
accounts for the difference between the group velocities 
of the vector components evolving freely (i.e. without 
coupling) , and 5 accounts for the gap between the energy 
dispersion curves of each component. These two later pa- 
rameters are of no fundamental importance, as we have 
demonstrated selftrapping also for the case An = 6 = 0, 
but should be taken into account to describe particular 
physical systems as BEC in optical lattices, for which 
each component belongs to a different energy band. In 
many aspects, Eqs. Illll'l) resemble to those describing 
the formation of solitary waves in light passing through 
a periodically twisted birrefringent fiber [Lq . In both 
cases, the vector components interact through linear and 
non-linear couplings, and velocity and energy mismatch 
between the vector components are taken into account. 
However, in contrast with Eqs. (11121 . the dispersion 
terms in |12| have the same signs and the non-linear 
potentials are attractive, i. e. dispersion balances the 
non-linearity in each equation. Up to our knowledge, 
this characteristic is shared by the generality of models 
for vector solitons. Therefore, the novelty of the system 
(11121) is the opposite sign in the dispersion terms for the 
two components. As discussed above, in this case the 
coupling is fundamental for the stabilization of the sec- 
ond component of the vector solitary wave. Fig. ^shows 
the possibility of self-trapping through the numerical in- 
tegration of Eqs. (111211 . The initial wavefunction is a 
Gaussian in one of the vector components and zero in 
the other, therefore a situation quite different from the 
two-component solitary wave in which it finally evolves. 
The parameters are taken as a = 0.775, 6 = 0.0II4, 
f3 = -0.014, 71 = 1.43, 72 = 0.56 and the norm of 
the initial wavefunction is 0.69, this particular values 
come from estimations from the Bose-Einstein experi- 
ment computed below. Figure ^a) and (b) correspond 
to the time evolution of the first and second component 
of the vector function, respectively. In both cases a tran- 
sient time is observed before the system reaches a stable 
configuration. During this transient, the population os- 
cillates between both vector components while radiation 
is emitted at both sides of the wavefunction. Matter ra- 
diation is associated to the transient and may have a 
fundamental role as a dissipation process that allows the 
system to converge to the stable solution 0, 0] . The 
final population enclosed in the solitary wave depends on 
this transient, which is more pronounced the greater the 
difference between the initial wavepacket and the final 
solitary wave is. 

Let us now focus to the experimental realization of this 
kind of vector solitary waves. To this end we will first 
demonstrate that Eqs. (111211 can describe the evolution 
of a Bose-Einstein condensate in an oscillating optical 
lattice. Afterwards, as a further proof, we will compute 
directly the time evolution of the condensate, as given by 
Gross-Pitaevskii equation. These ab initio computations 



(2) will demonstrate that, for conditions close to the exper- 
iment 0, vector solitary waves with mixed dispersion 
appear spontaneously. 

Following the philosophy of we will consider a 
Bose-Einstein condensate prepared near the edge of the 
Brillouin zone of an optical lattice. The condensate is 
tightly confined in the transversal direction by a trap- 
ping potential. Under these circumstances, the conden- 
sate dynamics can be considered one dimensional, pro- 
vided (/i/huj± — f)/2 << I 13], fi being the chemical 
potential and lu± the transversal frequency of the trap. 
The dynamical equation is then 



ih^(z,t) = {H 



H„i)$(z,t) 
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where H = p 2 /2m + U(z,t), H n i = g\n |£(z,t)| , with 
9id — Nha s /ma±, N being the number of atoms, a s 
the scattering length of the condensed element (in our 
case 87 Rb) and a± the transversal size of the wavefunc- 
tion (assumed to be in the lowest state of the radial har- 
monic potential). The potential U(z,t) corresponds to 
a time-dependent axial lattice. In our case this poten- 
tial includes an oscillatory motion of frequency u>, and 
amplitude £o- Therefore U(z,t) = Uo[z — £(i)] with 
C(^) — Cosin(wi), and Uq(<;) = Vq cos 2 (7r<r/ao) the sta- 
tionary lattice potential. 

Let us now define a new wavefunction </> according to 
the transformation 

<j>(z, t) = exp[i|c(*)]*(«. *) = ®(z + C(*)) (4) 



Eq. @ may now be written as 
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in which the time dependence is removed from the poten- 
tial term and translated to a linear coupling term. Note 
that this equation is formally identical to the NLSE of a 
ID system in interaction with a semiclassical electromag- 
netic field in dipole approximation, the effective vector 
potential being A(t) = (mc/q)((t)- This optical analogy 
is fruitful in understanding lattice oscillations as direct 
interband transitions, in the same fashion as electromag- 
netic waves produce optical transitions of electrons in 
crystals. This type of coherent transfer between bands 
has been experimentally characterized in 0] . 

Now we express the wavefunction 4>{z, t) as a superpo- 
sition of states of the first and second band. Using the 
Wannier basis [l7|. 



(z,t) 
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where fc n = ravo/h corresponds to the initial momentum 
of the wavefunction, located near the edge of the Bril- 
louin zone. We may now follow the standard procedure 
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for the effective mass approach in the form used, for in- 
stance, in ^5j. Substituting iJBJ in JSJ, and assuming 
negligible overlap between Wannier functions of neigh- 
bor sites (tight-binding approx.), we find the dynamical 
equations for the wavefunction envelopes 



in f f P (Z,t) = f p f p (t,t) 



2a. 
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gin ]T r mm f*(z,t)f n (tt)f m (t;,t) (7) 

with band labels {p, q,n,m} valued 1 or 2, £ defined as 
the continuous limit of i, and with 7r p . g = if dza*(z — 
£)(—ih)da q (z — £ — ao)/dz and T p , q , n ,m = J dza*(z — 
£)a*(z — t)a n {z — t)a m {z — £) constant parameters. The 
kinetic operator is defined as 



T p = e p (k ) + v g ,p C 



-in-T— 



d 2 



dS, ) 2m* p <9£ 2 
(8) 

with e p (ko) the energy of the p band evaluated 
at the point fco of the Brillouin zone and v 9lP = 
(l/h)(d/dk)e p \ ko and 1/m* = (l/h 2 )(d 2 /d 2 k)e p \ k , the 
group velocity and the effective mass. The linear 
coupling in Eq. describes the direct (optical) 
interband transitions induced by the oscillating lat- 
tice. Its derivation follows from the evaluation of 
the linear coupling in Eq. JSJ projected on the 
Wannier function exp (iko£)a p (z — £) which leads to 
sum of the type J2n,e> ex P i^o(£' ~ £)] fn(£' ', t) J dza*(z — 
£)pa n (z — £'). Assuming that the leading term in 
the sum corresponds to the coupling between near- 
est neighbors, and since fcoao ~ 7r, we can approxi- 
mate this sum to -i J2 n n P ,n [fn(£ + ao,t) - f n (l - a Q , t)] 
which, in the continuous limit, can be written as 

We shall assume that the lattice oscillation is tuned 
near the resonance of the interband transition: hu — 
e 2 (fc ) - m 2 w 9,2/ 2 - ei(fco) + m\v 2 g l /2 - HA, A being a 
small qua ntity. Following the standard approach to reso- 
nance |l8j . we define new rotating system wavefunctions 
g p according to 
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with Hk = (m*Vg t i + m2f 9j 2)/2. In the philosophy of the 
rotating- wave approximation |18| . we discard oscillating 
terms after substituting (|9I1U|) into Q, thus retaining 
only the slowly varying dynamics. Under these approxi- 
mations, Eq. (7J can be finally cast to the form of Eqs . 



Q]-|2J defining adimensional variables r] = £/ a o and r = 
i/ro, with To = y^mfrof|ao/?i, and defining the parame- 
ters as a — \l\m\lmX\, Ak = ao(m2Vg } 2—rnlv g ^i)/h, 5 = 
At , = -aoToAQWic lj2 k/K, and Uj = v / 2g 1D Ti :2 T( ) /hg j 
, and 7i = Ti^i/Ti^z and 72 = r 2 ,2 : 2,2/I\i,2,2- 

To check the validity of our approach, we have also 
performed ab initio integrations of the evolution of a 
Bose-Einstein condensate in an optical lattice, as given 
by Eq. J2J . In these computations, the transversal trap- 
ping is assumed to be constant in time and harmonic with 
uj± = 2ir x 50Hz. Initially a condensate of 400 atoms is 
prepared in a trap formed by this transversal potential 
and an axial potential with wm — 2 x 9Hz, (point (a) 
in figure |2J) ■ At the starting time, the axial potential 
is switched off, and the condensate is phase imprinted 
with a momentum hko = h0.9ir/ao with clq — 1.96/xto. 
In the experiment Q this phase imprinting is replaced 
by a displacement of the optical lattice with constant 
velocity. After the phase imprinting, an oscillating op- 
tical lattice potential is switched on adiabatically dur- 
ing 10ms to reach a maximum potential of Vq — 0.85E r 
(E r = h 2 Tr 2 /2mao). As the lattice period is ao, the 
wavepacket imprinted with momentum hko is now lo- 
cated near the edge of the first Brillouin zone (points 
(b) and (c) in figure |3J). The frequency of the lattice 
oscillation is u> = 2ir x 64Hz and the oscillation ampli- 
tude £0 = 0.35/im. With this procedure we obtain a 
wavepacket with similar characteristics that in the ex- 
periment but in an oscillating lattice. Figs|3fa) and 
(b) show the results of this ab initio integration of Eq. 
(J3J) in a similar fashion as in plots ^a) and (b) respec- 
tively. In Fig. |3 the vector components /1 and f 2 have 
been extracted from the time-dependent wavefunction by 
projecting in the Wannier base of the first and second 
bands. The mechanism of formation of the solitary wave 
follows the same steps as in the model eqs. QH2[1 : a 
transient characterized by the emission of radiation, and 
final state in which the two-component solitary wave is 
formed. Note that, in contrast with Fig. ^ no oscilla- 
tions between the vector components are observed. This 
is not a general case, as our computations show that the 
oscillating transient behavior can be recovered when the 
form of the initial wavefunction is chosen to be more sim- 
ilar to the solitary wave shape. We have found, in this 
case, that the frequency of the oscillations increases with 
the amplitude Co and as frequency to decreases (i.e. S in- 
creases). This suggests that these oscillations are of the 
Rabi type. The absence of them in Fig. [3] can be then 
explained by the more complex dynamics involved in the 
convergence of the system as the difference between the 
shape of the initial wavefunction and the final solitary 
wave is larger. For the case of scalar gap solitons, it is 
known that the wavepackets drifts in momentum space 
as a consequence of the radiation process during the tran- 
sient llU- Under these circumstances, the parameters of 
the interband transitions change and, therefore, the Rabi 
regime is severely distorted. 

In conclusion, we have proved the existence of a new 
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kind of vector solitary wave in which the dispersion sign 
is opposite for each component. Under these circum- 
stances, one of the vector elements is highly dispersive. 
Selftrapping is obtained through a linear coupling in the 
vector equations, which allows the highly dispersive com- 
ponent to be stabilized. We have demonstrated that this 



kind of vector solitary waves can be found in experimen- 
tally using Bose-Einstein condensates in oscillating opti- 
cal lattices. 
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FIG. 1: Time evolution of a two-component wavepacket as 
described by the coupled system 11121 . The wavepacket is 
initially a Gaussian in the lowest component (a) and empty 
in the upper component (b). After a transient time in which 
population oscillates between both vector components and 
radiation takes place, a two-component vector solitary wave 
is reached. The vertical scale in (b) has been magnified a 
factor two form that of (a). 



FIG. 2: Band structure of the optical lattice considered in the 
text. The condensate is initially located at momentum (a). 
After imprinting a momentum hko and the adiabatic switch 
on of the oscillating lattice, the wavefunction is found located 
at the edge of the Brillouin zone with components in the two 
lowest bands (b) and (c). 



FIG. 3: Time evolution of the condensate's wavepacket 
as computed from the Gross-Pitaevskii equation The 
wavepacket at t = corresponds to the lowest eigenstate of 
this equation for the trap defined in the text. At t = the con- 
densate is imprinted with a momentum hko. Subsequently an 
oscillating optical lattice is grown adiabatically. Plot (a) and 
(b) show the projection of the wavefunction onto the Wannier 
basis of the first and second band respectively. The vertical 
scale in (b) has been magnified a factor two form that of (a) 
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